Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



Constraints on the emission mechanisms of gamma-ray 
bursts 



on 

ON 
ON 

o 
Q 

(N 



(N 

ON 
ON 

Oh 
6 

C3 



Gabriele Ghisellini 1 , Annalisa Celotti 2 and Davide Lazzati 1 ' 3 , 

1 Osservatorio Astronomico di Brera, Via Bianchi 46, 1-23807 Merate (Lc), Italy 

2 SISSA, Via Beirut 2-4, 1-34014 Trieste, Italy 

3 Dipartimento di Fisica, Universitd degli Studi di Milano, Via Celoria 16, 1-20133 Milano, Italy 
E-mail: gabriele@merate.mi.astro.it, celotti9sissa.it, lazzatiamerate.mi.astro.it 



1 February 2008 



ABSTRACT 

If the emission of gamma-ray bursts were due to the synchrotron process in the stan- 
dard internal shock scenario, then the typical observed spectrum should have a slope 
F v oc v^ 1 / 2 , which strongly conflicts with the much harder spectra observed. This 
directly follows from the cooling time being much shorter than the dynamical time. 
Particle re-acceleration, deviations from equipartition, fastly changing magnetic fields 
and adiabatic losses are found to be inadequate to account for this discrepancy. We 
also find that in the internal shock scenario the relativistic inverse Compton scattering 
is always as important as the synchrotron process, and faces the same problems. This 
indicates that the burst emission is not produced by relativistic electrons emitting 
synchrotron and inverse Compton radiation. 



Key words: 

thermal 



gamma rays: bursts — X-rays: general — radiation mechanisms: non- 



1 INTRODUCTION 

Since the observational breakthrough by BeppoSAX (Costa 
, et al. 1997; van Paradijs et al. 1997) the physics of gamma- 
ray bursts (GRB) has started to be disclosed. The huge en- 
ergy and power releases required by their cosmological dis- 
tances support the fireball scenario (Cavallo & Rees 1978; 
Rees & Meszaros 1992; Meszaros & Rees 1993), whose evo- 
lution and behavior is (unfortunately) largely independent 
of their origin. 

We do not know yet in any detail how the GRB event is 
related to the afterglow emission, but in the most accepted 
picture of formation of and emission from internal/external 
shocks (Rees & Meszaros 1992; Rees & Meszaros 1994; Sari 
& Piran 1997), the former is due to collisions of pairs of 
relativistic shells (internal shocks), while the latter is gen- 
erated by the collisionless shocks produced by shells inter- 
acting with the interstellar medium (external shocks). The 
short spikes (t var ~10 ms) observed in the high energy light 
curves suggest that shell-shell collisions occur at distances 
R ~ 10 12 -10 13 cm from the central source, involving plasma 
moving with bulk Lorenz factor F > 100. The fireball starts 
to be decelerated by the interstellar medium further out, at 
a distance which depends on the density of this material. 

The main radiation mechanism assumed to be respon- 
sible for both the burst event and the afterglow is syn- 
chrotron (Rees & Meszaros 1994; Sari, Narayan & Piran 
1996; Sari & Piran 1997; Panaitescu & Meszaros 1998 - see 



however Thompson 1994; Liang 1997; Ghisellini & Celotti 
1999; Celotti & Ghisellini 1999; Stern 1999). This requires 
acceleration of electrons up to ultra-relativistic energies and 
the presence of a significant magnetic field. Evidence sup- 
porting that the afterglow emission is due to the synchrotron 
process include the power law decay in time of the afterglow 
flux (for reviews see Piran 1999; Meszaros 1999) and the re- 
cently detected linear polarization in GRB 990510 (Covino 
et al. 1999; Wijers et al. 1999), but the only piece of circum- 
stantial evidence in favor of a synchrotron origin of the burst 
radiation comes from the predicted frequency of the peak of 
the burst spectrum. Indeed, it is remarkable that the simple 
assumption of equipartition among protons, electrons and 
magnetic field energy densities leads - in the internal shock 
scenario (ISS) - to a typical emission frequency in agreement 
with observations. 

However, the very same ISS inevitably predicts very 
fast radiative cooling of the emitting particles. In this Let- 
ter we point out that this implies an emitted spectrum much 
steeper than observed. Although other authors have already 
pointed out that, in the presence of radiative losses, the 
predicted spectrum is steep (Cohen et al. 1997; Sari, Piran 
& Narayan 1998; Chiang 1999), here possible alternatives 
to avoid this conclusion, in the context of the ISS, are dis- 
cussed, but found inadequate to account for the discrepancy. 
In addition, we examine the role of the relativistic inverse 
Compton process in the ISS, which results to be as impor- 
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tant as the synchrotron one. Therefore in this scenario the 
high energy radiation would always be energetically signifi- 
cant, thus requiring a careful estimate of the importance of 
photon-photon collisions leading to electron-positron pair 
production. 



2 THE 'STANDARD' SYNCHROTRON 
SCENARIO 

Let us briefly summarize the main features of the ISS (for 
simplicity we will also refer to it as the 'standard' model). 
The emission of the burst, which originates from the conver- 
sion of bulk kinetic energy into random energy, has a dura- 
tion which is determined by the central engine. In particular, 
in order to generate intermittent and complex variability 
patterns, the engine has to produce/eject several shells of 
matter. If these propagate with different Lorentz factors, a 
faster shell will catch up with a slower one and in the inter- 
action a shock will develop, which is assumed to be respon- 
sible for the acceleration of electrons to ultra-relativistic 
energies. These would then loose energy by radiating syn- 
chrotron photons. 

2.1 Typical radii 

If two shells move with different Lorentz factors, F and aT 
(with a > 1), and are initially separated by R a , they interact 
at a distance Ri from the central engine, where 



Ri 



2a 2 



RoT 2 



(1) 



As a reference value in the following we will adopt a = 2, 
corresponding to Ri ~ (8/3)i? o r 2 . The initial temporal sep- 
aration of the two shells, ~ R a /c, also determines the dura- 
tion of the emission produced by a single shell-shell collision, 
as measured by an external observer. 

For a — 2 the final bulk Lorentz factor of both shells 
is ~ 1.4ir and if they have equal masses the difference be- 
tween the initial and final energy is ~6 per cent of their 
total kinetic energy. In the ISS, this energy is shared among 
magnetic field, protons and electrons. 

The other typical distance characterizing the fireball 
evolution is the transparency radius Rt, at which an ex- 
panding shell becomes optically thin to Thomson scattering. 
Assuming that each shell carries an energy E B — 10 50 i?a,50 
erg in bulk motion, R t ~ 6 X I0 12 (£: s ,5o/r2) 1/2 cm[]. 

The standard scenario requires Ri > R t , i.e. 
r ^ 350(-E Sl 5o/i?o,7) 1,/5 ! m order for the radiation produced 
to freely escape (see e.g. Lazzati, Ghisellini & Celotti 1999). 
In the following E s and L s stand for the kinetic energy and 
power of each shell, respectively, while L stands for the ob- 
served luminosity. 



Lorentz factor which corresponds to equipartition with the 
other forms of energy, i.e. 7 cq = (r' — l)n p m p / (n e m c ), where 
r' is the Lorentz factor of one shell in the rest frame of the 
other, and n p and n e axe the proton and lepton densities, re- 
spectively. In the ISS, these are assumed equal (i.e. electron- 
positron pairs do not significantly contribute to n c ). Devi- 
ations from the equipartition value are parametrized by a 
dimensionless coefficient e c = 7/7 cq . 

The out-flowing plasma is magnetized, and a typi- 
cal/indicative field value is estimated by assuming that ei- 
ther a significant fraction of the power is carried as Poynting 
flux or the field energy in the emitting region constitutes 
some fraction ee of the randomized energy. Both possibili- 
ties imply that at the distance where the shells interact the 
Poynting flux carries a power Lb = R 2 T 2 B' 2 c/2 = €b£ s , 
where L B = 4irR 2 r 2 n p m p c i is the kinetic power carried by 
a single shell. This corresponds to 



B oq = \8neBn p m p c \ = I — - — I 



2e B L s \ 1 / 2 J_ 
TR 



(2) 



From the above estimates it follows that the typical observed 
synchrotron frequency is f pca k = 2e/(37rm c c)7 2 B'r/(l + z), 
which at Ri gives 



\2 1/2,1/2 



Ri,i 3 (l + z) 



MeV. 



(3) 



As mentioned in the Introduction, the success of the stan- 
dard scenario in (simply) predicting a typical observed fre- 
quency in remarkable agreement with observations is prob- 
ably the strongest piece of evidence pointing towards the 
synchrotron process as responsible for the burst emission. 

Note that the 'equipartition coefficients', e B and e e , 
must be close to unity Jfor the observed value of f pca k to be 
recovered (note that (r' — 1) is of order unity in the ISS). In 
turn this also implies/requires that electron-positron pairs 
cannot significantly contribute to the lepton density. 

The predicted synchrotron spectrum, produced by a 
quasi mono-energetic particle distribution, has a flux den- 
sity F v oc i/ 1 / 3 up to the cutoff frequency f pca k. The average 
observed spectra in the hard X-ray band, which we stress 
are typically derived from ~ 1 s integrated fluxes, are not 
inconsistent with this shape. Nevertheless exceptions exist, 
including spectra much flatter than v x l' A , which have already 
cast some doubts on the synchrotron scenario (Preece et al. 
1998; Lloyd & Petrosian 1999). 

In the following we point out that, just because the 
integration and the dynamical timescales are much longer 
than the particle cooling timescales, the expected synchrotron 
spectrum in the entire X- and soft "/-ray band should have 
a slope F v oc u^ 1 ^ 2 . This dramatically exacerbates the dis- 
crepancy between the predictions of the standard scenario 
and observations. 



2.2 The typical synchrotron frequency 

During the shell interaction electrons are instantaneously 
accelerated in a collisionless shock, and reach a random 

* Here and in the following we parametrize a quantity Q as Q = 
10® Qx and adopt cgs units. Primed quantities are evaluated in 
the comoving frame. 



T The only estimates of these parameters are inferred from the 
adoption of this model for the interpretation of the afterglow emis- 
sion. This suggests values of ee and e e substantially smaller than 
unity. 
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3 RADIATIVE COOLING TIME AND TIME 
INTEGRATED SPECTRUM 

Consider the radiative cooling timescale (in the observer 
frame) of typical particles emitting synchrotron (and self- 
Compton) radiation within the frame of the ISS: 



, _ 7 1+z _ 

£cool — . - 

7 r 

= 1.14 x 10" 



67rm e c(l + z) 

<j T B 2 r 7 (i + Uy/Ub) 
- i) 3 r 2 

<,V(1 + Ur/U B )(1 + Z) 



(4) 



where U T and Ub represent the radiation and magnetic en- 
ergy densities, respectively. As already mentioned, the short- 
est integration times are of the order of 1 s: this implies 
that the observed spectrum is produced by a cooling parti- 
cle distribution. Note also that the cooling timescale is much 
shorter than the dynamical one, resulting in a relatively effi- 
cient radiative dissipation: in this situation adiabatic energy 
losses are therefore negligible (Cohen et al. 1997; see below). 
In particular, after one dynamical time t<j = 10™ 2 id,-2 s, the 
cooling electrons emit at the (cooling) observed frequency 
^cooi ~ 2.2 x 10 14 (1 + 2)td 2 _ 2 r" 1 B~ 3 (1 + Ur/U B )- 2 Hz, in- 
dependent of r'. Invoking smaller values of the magnetic 
field to slow down the cooling and obtain v coo \ of the order 
of few hundreds keV does not help, since in this case the 
self-Compton emission dominates over the synchrotron one 
(see below). 

Since tcooi oc I/7, in order to conserve the particle num- 
ber, the instantaneous cooling distribution has to satisfy 
N(-y,t) oc I/7. When integrated over time, the contribution 
from particles with different Lorentz factors is 'weighted' 
by their cooling timescale oc I/7. Therefore the predicted 
(time integrated) flux spectrum between v coo \ and 
(e.g. Piran 1999) 



F v oc t coo i N(j) 7 — oc — 7 v 
dv 7 7 



1/2 -1/2 



(5) 



extending from ~ hv CO o\ ~eV to /i^ pC ak ~MeV energies. We 
thus conclude that, within the assumptions of the ISS, a 
major problem arises in interpreting the observed spectra 
as synchrotron radiation. Let us consider in turn alterna- 
tive hypotheses, within the same general frame, which might 
ameliorate this difficulty. 



3.1 Deviations from equipartition? 

If one maintains the requirement of observing synchrotron 
photons at ~MeV energies, the radiative cooling timescale 
is almost independent of eb, i.e. of the assumed value of the 
magnetic field (the only dependence being through the ratio 
U t /Ub; eq. 4). Furthermore - again from eq (4) - t coo i ~ i d 
requires a value of e c close to ~ 40, thus violating energy 
conservation (as the electrons would have more energy than 
the available one). 

On the other hand, if the condition that ^ pea k ~1 MeV 
is produced by synchrotron is relaxed, the magnetic field 
intensity can be smaller than the equipartition value with a 
consequently longer synchrotron cooling timescale. However, 
as the radiation energy density has to be of the order of U[ ~ 
L s /(A-kR 2 T 2 c) ~ f]3,cq to account for the observed fluxes, 
the inverse Compton cooling would be in any case extremely 
efficient, leading again to short cooling timescales. Therefore 



even if the observed radiation is produced by self-Compton 
emission of relativistic particles, we face the same problem 
of cooling timescales being so short that the spectrum would 
be steep, as discussed below. 

3.2 Particle re-acceleration? 

A further possibility to escape the above conclusion is to as- 
sume that particles are continuously re-heated, thus avoid- 
ing the formation of a cooled particle distribution. However 
in the standard ISS new particles are continuously swept by 
the shock and are all accelerated to the equipartition energy. 
This is a crucial assumption in order to produce a typical 
observed peak frequency around a few hundred keV. It is 
thus not possible - in this scenario - to continuously re- 
accelerate the very same particles, as the energy required 
would exceed the available one. 

Alternatively, relaxing the requirement of the standard 
scenario, one can envisage a situation in which only 'selected' 
particles are steadily accelerated for the entire duration of 
the shell-shell interaction. In this case an extreme fine (and 
unlikely) tuning is required: in fact, to be consistent with the 
total energetics, the selected particles have to: i) be fixed in 
number (only a fraction ~ t' cool /(AR' /c) of the total number 
of particles can be accelerated); ii) be always the same; iii) 
achieve 7 ~ 7 eq even in the absence of an equipartition 
argument. 

It would be also plausible to assume that the emission 
is produced by a power-law distribution of electrons result- 
ing from continuous acceleration and cooling. Indeed, for an 
energy distribution oc j~ p (with p > 0) only a minority of 
particles attain the maximum energy. But - besides having 
to keep all the particles accelerated for the entire duration of 
the shell-shell interaction - the relative number of the most 
energetic particles requires p > 2, leading to a spectrum 
even steeper than u^ 1 ^ 2 . 

We therefore conclude that re-acceleration does not 
avoid the spectral discrepancy, even when relaxing some of 
the key assumptions of the standard scenario. 



3.3 Strongly varying magnetic field? 

Let us consider the case in which the magnetic field attains 
a value close to the equipartition one only in a very lim- 
ited region (e.g. near the shock front), while is weaker else- 
where. In this situation the synchrotron cooling is mostly 
effective within this region only and the particles may not 
have time to significantly cool. Therefore in principle a syn- 
chrotron spectrum oc v 1 '' 3 might be produced. This requires 
that particles loose much less than half of their energy in 
the radiative zone, since even a reduction of a factor two 
of their Lorentz factor would imply that a spectrum u^ 1 ^ 2 
is produced in a range spanning a factor four in frequency 
(this may correspond to the entire BATSE energy range). 
Therefore it would be required that: 

i) the synchrotron process in the most radiative region must 
be inefficient, since it has to reduce the electron energy at 
most by a small fraction; 

ii) away from this zone, particles continue to rapidly cool 
by self-Compton and - at a reduced rate - by synchrotron 
emission. The inverse Compton process then becomes the 
dominant cooling mechanism. 
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iii) since the cooling is very rapid anyway, the self-Compton 
emission itself would produce a time integrated (over id) 
steep spectrum. 

We conclude that the net effect of having a strong mag- 
netic field confined in a limited region is to decrease the total 
synchrotron power in favor of the self-Compton one, whose 
spectrum would in any case be steep (see below). 

3.4 Adiabatic losses? 

Suppose that particles are accelerated in compact regions 
that rapidly expand because of internal pressure. Adia- 
batic losses dominate particle cooling as soon as the particle 
Lorentz factor decreases below some critical 7 a( j, thus gen- 
erating a spectrum oc u 1 ^ 3 below the synchrotron frequency 
^ a d (corresponding to 7 a a), an d steeper above (Cohen et al. 
1997). However this possibility faces two severe problems, 
both related to the overall efficiency, being required that: 

i) each electron loses only a small fraction of its energy ra- 
diatively (i.e. 7ad/7 must be greater than ~ 1/2); 

ii) the emitting regions are very compact, for adiabatic losses 
to be significant. This implies that the transformation of 
bulk into random energy does not occur in a shell subtend- 
ing the entire ejection solid angle. Photon and electron den- 
sities have then to be higher to account for the observed 
luminosity, thus enhancing the inverse Compton process. 



4 IMPORTANCE OF THE RELATIVISTIC 
INVERSE COMPTON PROCESS 

As long as the scattering optical depth tt of the electron in 
the emitting region is smaller than unity, the importance of 
the relativistic inverse Compton process with respect to the 
synchrotron one is measured by the (relativistic) Compton 
parameter y' = tt7 2 /3 2 = <tt7 2 (3 2 n' e ct' cool . The width of the 
region corresponds to a cooling length, ct' cool , as assumed 
within the standard ISS. 

Since the magnetic field intensity is related to the pro- 
ton density n' p we obtain 

/ 3 e c n c r 1 , . 

y ~ I ^ < i + m/u^ ' 1 ' 

This implies that the inverse self-Compton power is of 
the same order of the synchrotron one [J], and is emitted at 
a typical observed energy (for the first order) 



7 2 /w a 



13 



4(r' 



i) 4 4 /2 ^ 5 2 2 



fli,is(l + z) 



TeV. 



(7) 



Two points are worth being stressed. First, the strong de- 
pendence of u c on the equipartition parameter e e - Further- 
more, although in the comoving frame this typical Compton 
frequency is a factor T lower, it can still largely exceed the 
pair production threshold (see below). 

It has been mentioned and implicitly assumed above 
that even in the hypothesis that the hard X-ray burst ra- 
diation is due to self-Compton emission, the argument of 
the fast cooling producing a steep spectrum applies. Let us 

T This is true as long as the scattering process is in the Thomson 



consider this possibility more closely, and in particular the 
first order inverse Compton spectrum. 

Although in this case the typical electron energies re- 
quired are smaller, the cooling timescales are still much 
shorter than the dynamical time: in fact, to produce 
~MeV photons by the first order Compton scattering, 7 ~ 
83[^Mcv(l+^)/(r2-B4)] 1 ^ 4 with a corresponding cooling time 



t cool = 1.4 x 1(T 5 -^-^ 4 !: — s 



N3/4 



(8) 



L 50 (l + UB/U r )^ v 

Here L = 10 50 L5o erg s _1 is the (observed) luminosity pro- 
duced by a single shell. 

Furthermore, by following the same arguments leading 
to eq. (5), the predicted time-integrated spectrum results 
F v oc i/ -3 / 4 , i.e. even steeper than v~ x l' 1 . 

A further difficulty of interpreting the burst emission as 
first order scattering is that if the inverse Compton power 
exceeds the synchrotron one by a certain factor, then each 
higher Compton order will dominate over the previous one 
by the same amount, until the typical emitted frequency 
reaches the electron energy. Only a small fraction of the 
radiated power would therefore be observed (in the hard 
X-ray band). 



5 PAIR PRODUCTION 

The above results indicate that the time integrated spectrum 
predicted by the standard scenario is steeper than observed. 
Furthermore, the power emitted through the self-Compton 
process should be comparable to - if not more than - the 
synchrotron one, and emitted at energies exceeding the pair 
production threshold. It is thus compelling to estimate the 
importance of photon-photon collisions producing electron- 
positron pairs. 

Setting x = hu/(m c c 2 ), the energy threshold for pho- 
tons of energy x is xt — 2/[x(l — cos#)], where 9 is angle 
between the two photon directions. Also, the photon-photon 
collision rate is proportional to (1 — cos 9). 

The result of the integration of the photon-photon cross 
section over the energy of the target photons can be well 
approximated by (<jt /5)xttIj(xt) (Svensson 1987), where 
n-y(xT) is the number density of photons of energy xt, which 
is related to the observed luminosity L(xt) by XTn-,(xT) = 
L(x T )/(47rm c c 3 i? 2 ). 

The optical depth for pair production in the observer 
frame can be then expressed as 



_ °t L(x T )(l - cos 9) AR 
Tll[X) ~ 20^ 7W 3 W 



(9) 



regime, i.e. 7 < 800B. 



-1/3 



for the first order Compton scattering. 



AR may represent the width of the emitting shell or, alter- 
natively, the typical scale over which the emitted photons 
might interact (i.e. AT? ~ R), depending whether we are 
interested in the pair production within the shell or also 
outside it. 

Since the source is moving relativistically, all photons 
appear to be emitted quasi-radially, and interact with a typ- 
ical angle sin 9 ~ 1/T (corresponding to cos 9^/3), and thus 
(1 - cos0) ~ 1/r 2 . If the typical size of the fireball is esti- 
mated by time variability, 7? ~ ct vs , r F 2 , we have: 

L(x T ) AR 



I \ — °" T 



20tt i var m c c 4 r 4 R 



(10) 
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In this form r 77 (a;) can be estimated even without detailed 
spectral information. If the observed spectrum is a power 
law L(x) oc x~ a , with a < 1 up to a maximum energy 
i m „, the observed luminosity at threshold is related to the 
total luminosity L by L(xt) = [(1-a) /2 a ](L/x^£)(x/T 2 ) a , 
giving 



(x/2) a (1 -a)l AR 
T 4 + 2 - 207rai=S R ' 



(11) 



where the compactness £ = aT L / (t va , T m e c 4 ) has been intro- 
duced. Note that the optical depth increases with photon 
energy x. 

For illustration, consider a burst with V — 10 2 lasting 
tvar = 10 ms. 1 GeV photons (x = 2000) mostly interact 
with target photons of energies xt = 10 rjj, i.e. ~ 5 MeV. 
Assume that the observed luminosity at these energies is 
L(x T ) = 10 50 erg s" 1 . From eq. (10) we have 



r 77 (a; = 2000) ~ 1.4 X 10 : 



L 50 (x T ) 

1 2 



AR 
R ' 



(12) 



Since the optical depth is so large, all the high (~GeV) en- 
ergy emission can be easily absorbed, unless T > 10 3 . 

For r < 10 3 relativistic pairs can then be copiously pro- 
duced. They will immediately cool radiatively, initiating a 
pair cascade, strongly affecting the primary spectrum, and 
possibly even the dynamics. (Equilibrium) spectra produced 
by pair cascades have been extensively studied in the past 
in the context of nuclear AGN emission. The general out- 
come is that pairs act as reprocessors of the high energy 
emission, which is absorbed and ultimately reprocessed into 
lower energy radiation. This corresponds to a steepening of 
the spectrum, thus exacerbating the discrepancy with the 
observed bursts. 

Furthermore, the lepton density n e may become sub- 
stantially larger than the proton one, n p , and thus the 
equipartition energy 7 cq smaller than m p /m s . 



6 DISCUSSION 

The main point stressed in this paper is the inadequacy of 
the synchrotron and inverse Compton emissions from ultra- 
relativistic electrons to account for the observed burst spec- 
tra - at least within the scenario invoking internal shocks 
for the dissipation of the fireball bulk kinetic energy. 

This is independent of detailed assumptions and di- 
rectly follows from the extremely short cooling timescales 
- compared to the dynamical time which is in turn much 
smaller than the integration one - required by the ISS, which 
lead to a steep emitted spectrum. No alternative hypothe- 
sis, which could alleviate this spectral discrepancy, has been 
found. Furthermore it is stressed that within the ISS sce- 
nario electron-positron pairs would be naturally and copi- 
ously produced, contrary to the basic model assumptions. 

The situation is somewhat paradoxically: the emission 
mechanism at the origin of the burst radiation must be very 
efficient, and the synchrotron and inverse Compton mecha- 
nisms by relativistic particles are indeed very efficient radi- 
ation processes, but just because of the very rapid cooling 
their predicted spectrum is too steep. 

One is therefore forced to look for alternatives. In the 
dense photon environment of the internal shock scenario a 



highly efficient viable alternative radiation mechanism may 
be Comptonization by a quasi-thermal particle distribution, 
as proposed by Ghisellini & Celotti (1999) (see also Thomp- 
son 1994; Liang 1997; Stern 1999; Liang et al. 1999). In this 
model, the conversion of bulk kinetic into random energy 
may still be due to shell-shell collisions, but the typical en- 
ergy of the radiating particle is sub-relativistic, being fixed 
by the balance between the acceleration and the cooling pro- 
cesses. There is still equipartition between magnetic field, 
leptons and protons energies, but in a time integrated sense: 
all leptons are accelerated up to small energies, but for the 
entire duration of the shell-shell collision. Electron-positron 
pairs can be produced, and may even be the key ingredient 
to lock the particle energies in the observed range. These 
sub- or mildly relativistic particles would then emit self- 
absorbed cyclo-synchrotron photons and a Comptonization 
spectrum with a typical slope F v oc v° plus a Wien peak 
located where photon and particle energies are equal. The 
predicted spectral and temporal evolutions from this model 
are under investigation. 
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